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Theoretical analysis is given of nonclassicality and decoherence of the field states generated by 
adding any number of photons to the squeezed thermal state (STS). Based on the fact that the 
squeezed number state can be considered as a single-variable Hermite polynomial excited state, the 
compact expression of the normalization factor is derived, a Legendre polynomial. The nonclas- 
sicality is investigated by exploring the sub-Poissonian and negative Wigner function (WF). The 
results show that the WF of single photon-added STS (PASTS) always has negative values at the 
phase space center. The decoherence effect on PASTS is examined by the analytical expression of 
WF. It is found that a longer threshold value of decay time is included in single PASTS than in 
single-photon subtraction STS. 
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I. INTRODUCTION 

Generation and manipulation of non-classical light 
field has been a topic of great interest in quantum op- 
tics and quantum information science [Ij]. Many experi- 
mental schemes have been proposed to generate nonclas- 
sical states of optical field. Among them, subtracting 
photons from and/or adding photons to quantum states 
have been paid much attention because these fields ex- 
hibit an abundant of nonclassical properties and may give 
access to a complete engineering of quantum states and to 
fundamental quantum phenomena [3-[l0j. For example, 
quantum-to-classical transition has been realized exper- 
imentally through single-photon-added coherent states 
of light. These states allow one to witness the gradual 
change from the spontaneous to the stimulated regimes 
of light emission [J] . For m- photon-added coherent state 
in the dissipative channel, the nonclassical properties are 
studied theoretically [11] by deriving the analytical ex- 
pression of the Wigner function (WF), which turns out 
to be a Laguerre-Gaussian function. As another exam- 
ple, photon addition and subtraction experimentally have 
been employed to probe quantum commutation rules by 
Parigi et al. In fact, they have implemented simple al- 
ternated sequences of photon creation (addition) and an- 
nihilation (subtraction) on a thermal field and observed 
the noncommutativity of the creation and annihilation 
operators [6]. In addition, photon subtraction/addition 
can be applied to improve entanglement between Gaus- 
sian states jH, nil , loophole- free tests of Bell's inequality 
[3, [Tsj , and quantum computing [l6| . 

On the other hand, it is interesting to notice that sub- 
tracting or adding one photon from/to pure squeezed vac- 
uum can generate the same output state, i.e., squeezed 
single-photon state fl7|. Actually, the photon addition 
is able to generate a nonclassical state (e.g coherent and 
thermal states), which is quite different frorn. photon sub- 
traction only from a nonclassical state |18l - [20| . In addi- 



tion, the resulting states obtained by successive photon 
subtractions or additions are different from each other. 
For instance, successive two-photon additions [a^"^] and 
successive two-photon subtractions [a^] will result in the 
same state produced by using subtraction- addition (a^a) 
and addition-subtraction (aa^), respectively. In Ref . j2]| . 
two photon-subtracted squeezed vacuum is used to gen- 
erate the squeezed superposition of coherent states with 
high fidelities and large amplitudes. 

In general, different non-Gaussian operators (e.g sub- 
tracting and adding photon) will suffer different effects 
from the surroundings, thus it is important to know 
which operator is more robust compared to the other 
under an identical initial quantum state when the envi- 
ronment is taken into account. Very recently, the robust- 
ness of several superposition states is studied by using the 
linear entropy under a thermal environment [22| . In this 
paper, we shall introduce a kind of nonclassical state — 
photon-addition squeezed thermal state (PASTS), gener- 
ated by adding photon to squeezed thermal state (STS) 
which can be considered as a generalized Gaussian state. 
Then we shall investigate the nonclassical properties and 
decoherence of single-mode any number PASTS under 
the influence of thermal environment. 

This paper is organized as follows. In Sec. II we intro- 
duce the single-mode PASTS. By converting the PASTS 
to an Hermite polynomial excitation squeezed vacuum 
state, we derive a compact expression for the normal- 
ization factor of PASTS, which is an m-order Legendre 
polynomial of squeezing parameter A and mean number 
Tic of thermal state, where m is the number of added pho- 
tons. In Sec III, we discuss the nonclassical properties of 
the PASTS in terms of sub-Poissonian statistics and the 
negativity of its WF. We find the negative region of WF 
in phase space and there is an upper bound value of A 
for this state to exhibit sub-Poissonian statistics which 
increases as m increases. Then, in Sec. IV we derive the 
explicitly analytical expression of time evolution of WF of 
the arbitrary PASTS in the thermal channel and discuss 
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the loss of nonclassicality in reference of the negativity of 
WF. The threshold value of decay time corresponding to 
the transition of the WF from partial negative to com- 
pletely positive definite is obtained at the center of the 
phase space, which is independent of parameters A and 
He. It shown that the WF for single PASTS (SPASTS) 
has always negative value for all parameters A and ric if 
the decay time nt < ^ \n[{2Af + 2) / {2Af + 1)] (see Eq.^ 
below), where Af denotes the average thermal photon 
number in the environment with dissipative coefficient 
K. Comparing to the case of single photon-subtraction 
STS (SPSSTS), the decoherence time of SPASTS is 
longer. In this sense, the photon-addition non-Gaussian 
states present more robust contrast to decoherence than 
photon-subtraction ones. The reason may be that the 
amount of non-Gaussianity for SPASTS is larger than 
that for SPSSTS as presented in Sec. V. Conclusions are 
involved in the last section. 
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where \q) is the eigenstate of Q = (a + a^)/V2, Q\q) = 
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Thus, using Eq.(IS|) and the overlap relation 
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where Hn {q) is the single- variable Hermite polynomial 
then 5*1 \n) can be expressed as 



II. PHOTON-ADDITION SQUEEZED 
THERMAL STATE (PASTS) 



The rn-photon-added scheme, denoted by the mapping 
p — >■ a^"^pa™, was first proposed by Agarwal and Tar a 
[is"]. Here, we introduce the PASTS. Theoretically, the 
PASTS can be obtained by repeatedly operating the pho- 
ton creation operator on a STS, so its density operator 
is given by 



C— 1 im ot o ^ 



(1) 



where m is the added photon number (a non-negative 
integer), is the normalization constant to be de- 



termined, and Si = exp[A( 



,2 _ «t2 



)/2] is the single- 



mode squeezing operator with A being squeezing param- 
eter [H, 111] . pth is a single field mode with frequency uj 
in a thermal equilibrium state corresponding to absolute 
temperature T, whose the density operator is [25?| 
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(: : denoting antinormally ordering) which implies that 
the density operator pth can be expanded as 
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where ric = [exp(a;/(fcT)) — being the average pho- 
ton number of the thermal state pth and ks being Boltz- 
mann's constant. Eq.Q is useful for later calculation. 
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where we have set sechA = 2p/{p^ + 1), tanh A = 
[p^ — l)l{p^ + 1), and we have used 5*110) = 
sech^/^Aexp[— a^^/2tanhA] |0) as well as the generating 
function of iJ„ (q) [13]: 



Hn (q) 



Qn 
J' 



- exp {2qT - 



(8) 



Eq.(I71) indicates that the single-mode squeezed num- 
ber state Si \n) is actually a Hermite polynomial ex- 
cited squeezed vacuum state [2$]. Obviously, when n — 
0, Hq (q) — I, Eq.([7]) just reduces to single- mode squeezed 
vacuum. While for n = 1,2, noting Hi{q) ~ 2q and 
H2 (q) = 4g2 - 2, Eq.0 become 

S-ill) ^a^sechA ^i |0) , 

Si |2) = ^ {a^'^sech^X + tanh A) Si |0) , (9) 
v2 

respectively. It is interesting to notice that the single 
photon-added squeezed vacuum (PASV) is equal to the 
squeezed number state 5*1 |1), and the two PASV can 
be considered as a superposition of the squeezed number 
state 5*1 1 2) and the squeezed vacuum. 



A. Squeezed number state as a Hermite 
polynomial excited state 

Recalling that the single-mode squeezed operator 5*1 
has its natural expression in the coordinate representa- 



B. Normalization of PASTS 

To fully describe a quantum state, its normalization is 
usually necessary. Next, we shall employ the fact ([7]) to 
realize our aim. First, let us derive the normally ordering 
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form of STS ps = S\pthSi, which is convenient for further 
calculation of normalization. 

Using Eqs.Q and dT]), we can rewrite the STS ps as 
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where 5*1 (—A) = Si (A) and the vacuum projector 
|0) (0| =: exp [— a'^a] : is used. Further using the two- 
linear generating function of Hermite polynomial [1^ , 



E (^) (y) 



n=0 



1 



exp 



1 -i2 



(11) 



we can directly obtain the normally ordering form of STS, 
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By introducing a = {Q + iP) /^/2 and = {Q ~ iP) /^/2, 
Eg. dT^ can be put into another form 
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where T1T2 ~ \/A^ and 
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Eq. ([T2|) or (|T4|) is a compact expression of the STS, which 
is just a Gaussian distribution within normal order for 
operators Q and P [30j . 

Next, we shall derive the normalization factor for 
PASTS. Employing Eq.dn]), the PASTS reads as 
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Thus the normalization factor Ca.m is (1 = trpad) 
— f ^ up™ e-(i-B)l"l'+f ("•'+"') 
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where we have used the completeness relation of coherent 
state, and [(1 - Bf' - C^]"^ = rfrl = A, as weU as the 
integration formula [3l| 



j ^ exp (C |z|' + + 77^* + fz^ + gz*^') 



exp 



(18) 



whose convergent condition is Re(C ± / ± g) < 0, Rc(C2 — 
4/.g)/(C±/±ff)<0. 

Recal ling the newly found formula of Legendre poly- 
nomial [32I [33|. i.e.. 
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2/2 
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and noticing .t^ — 1 = AC^, together with a; = 
VA {I - B) = [A - ric [ric + 1)] /\/~A, we have 



_ (AC)"' 



-{l~B)st-s^ -e 



(20) 



which indicates that Ca,m is also just related to Legendre 
polynomial, and 



= cosh 2A + cosh A. 



(21) 



It is noted that, for the case of no-photon- addition with 
TO = 0, Ca.o = 1 as expected. Under the case of m- 
photon-addition thermal state (no squeezing) with B = 
ric + l, A = {ric + 1)^ , and P™ (1) = 1, then Ca,,„ = 
to! (ric + 1)™ • The same result as Eq.(32) can be found 
in Ref.lsl. 



III. NONCLASSICAL PROPERTIES OF PASTS 

In this section, we shall discuss the nonclassical prop- 
erties of PASTS in terms of sub-Poissonian statistics and 
the negativity of its WF. 
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FIG. 1: (Color online) The Q-parameter as the function of 
squeezing parameterr for different m = 0, 1, 2, 3, 4, 19, 20 with 
a small ric value. 



A. Sub-Poissonian nature of PASTS 

The nonclassicality of the PASTS can be analyzed by 
studying its sub-Poissonian distribution. Using Eq. ipO)) 
we can directly calculate: 



(at a) - 



C'a.i 



Ca,' 
Ca,m+2 



-1, 



(22) 
(23) 



Thus the Mandel's Q-parameter [35l| can be obtained by 
substituting Eqs.dH]) into Q = (a^^a^) / (a^a) - (a^a) , 



Q = 



Ca. 



m+2 



Ca-, 



Ca.; 



(24) 

The negativity of the Mandel's Q-parameter refers to 
sub-Poissonian statistics of the state. In order to see 
clearly the variation of Q-parameter with A and ric, we 
show the plots of Q-parameter in Fig.l, from which 
one can clearly see that, for a given small ric value, Q- 
parameter becomes negative (m ^ 0) when A is less than 
a certain threshold value which increases as m increases; 
while for to = or a large rtc, Q is always positive. This 
implies that the nonclassicality is enhanced by adding 
photon to squeezed state. Here, we should emphasize 
that the WF has negative region for all A and ric, and 
thus the PASTS is nonclassical. 



B. Photon-number distribution (PND) of the 
PASTS 



The photon-number distribution (PND) is a key char- 
acteristic of every optical field. For this purpose, we first 
calculate the PND of STS, then the PND of PASTS can 
be directly obtain. The PND, i.e., the probability of 
finding n photons in a quantum state described by the 
density operator p, is V{n) = {n\p\n) . So the PND of 



the STS is 



V{n) = {n\ S\pthSi \n) 



(25) 



Using the fact in ([7]) and the P-representation of pth (Ellj 
Ea. (p5)) can be directly written as 
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In a similar way to deriving Eg. (1^01) . using Eq. p^ we 
have 



where 
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(28) 



Eq.(I27| shows that the PND of STS is the Legendre 
polynomial of Bj^fD. In particular, when \ — A — 
[n^ + Xf and Bl^fD = 1,D = nl/{nc + l)'^, then Eq.^ 
becomes V{n) = n"/ (ric + l)", corresponding to the PND 
of thermal state [3J|. In fact, we can also check Eq.([57]) 
using the normalization condition. Note that the Legen- 
dre polynomial can also be defined as the coefficients in 
a Taylor series expansion [36j 



1 



= ^F„(:r)r, 



(29) 



n=0 



thus E„=o^(") = 1/V^(1 -2B + D) = 1 as expected. 

Next, we turn to present th e PND of PA STS. From 
Ea. (l?7|) and noting a^™ \n) = y/ (m + n)\/n\ \m + n) and 
a™ |n) = y^n\/(n — m)l |n — to), it then directly follows 



V2{n)^C-l (n|at"p,a™|n) 



ilC7lD(-' 



i)/2 



{n — my.^/A 



-Pn- 



(b/Vd) 



(30) 



Eq.dSni) is the PND of PASTS, a Legendre polynomial 
with a condition n ^ m which implies that the photon- 
number (n) involved in PASTS is always no-less than the 
photon- number (m) operated on the STS, and there is no 
photon distribution when n < to). For some other non- 
Gaussian states, such as a'^"'a"^ p^a'^'^a" , a™a'^" psa"'a'^™ , 
and a"^psa^"^, their PNDs can also be directly obtained 
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FIG. 2: (Color online) Photon-number distributions of 
PASTS with n=l for A=0.3, m=0 (blue bar); A=0.3, m=l 
(red bar), A=0.3, m=5 (yellow bar), and A=0.8, m=l (green 
bar). 



by using In Fig. 2, the PND is shown for dif- 

ferent values (A, ric) and m. By adding photons, we have 
been able to move the peak from zero photons to nonzero 
photons (see blue and red bar in Fig. 2). The position of 
peak depends on how many photons are created and how 
much the state is squeezed initially. The probability of 
PND becomes smaller with the increasement of squeezing 
parameter (see red and green bar in Fig. 2). 



IV. WIGNER FUNCTION OF PASTS 



Next, the normally ordering form Eq.(|T2|) is applied 
to deduce the WF of PASTS. The partial negativity of 
WF is indeed a good indication of the highly nonclassical 
character of the state. Therefore it is worth of obtaining 
the WF for any states. The WF W{a,a*) associated 
with a quantum state p can be derived as follows [23| : 



W{a,a*) = e 



d2/3 



-/3|p|/3) e 



2(a/9*-a*/ 



(31) 



where |/3) = cxp(— \/3\^ /2 + /3a') |0) is the coherent state. 

Substituting Ea. (|T6|) into Eg. (1311) . we can finally obtain 
the WF of PASTS (see Appendix A), 

W (a, a*) ^F^ (a, a*) Wo (a, a*), (32) 

where Wq (a, a*) is the WF of STS, 
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and 
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(m!)^C-\sinh'"2A 
22™ {2n^ ^- 1)" 



^ ^![(m-0!] sinh'2A ' ' 



(34) 



where 7 = \q* sinh 2A — 2Q!(cosh^ A -I- 
nc)]/{i[(2nc + l)sinh2A]i/2}^ Eq.(|32l) is the analytical 
expression of WF for PASTS, related to single- variable 
Hermite polynomials. In particular, when m = 0, 
Fo (a,a*) = 1, Eq.dMl) becomes W {a,a*) = Wo («,«*); 
while for A = 0, note Cam — 'ti! (n^ + 1)™, Wo (a, a*) — 
e-2|«lV(2"c+i)/[^(2n^ + l)] and i^™(a,a*) 



(-l)'"/(2n, + l)"L, 
Eq.(I321) reduces to 



[4 (ne + l)|ar/(2nc + l)], 



W (a, a* 



i-iy 



e 2fi+i 



7r(2nc + I)'' 



4K + 1) , ,2 
'a 
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1 



(35) 

which corresponds to the WF of m-photon added thermal 
state [131, and can be checked directly from Eq.(A3). In 
addition, for m = 1, [single-photon-added squeezed ther- 
mal state (SPASTS)], Cai ^ B the special WF of 
SPASTS is 



Wi{a,a*) = Fi{a,a*)Wo{a,a*), 



where 
Fi{a,a*) 



sinh 2 A 

(2nc + 1)B 



cosh A 



I7l 



sinh 2 A 



(36) 



(37) 



Noting B > Q, thus from Eg. (1571) one can see that when 
the factor Fi (a, a*) < 0, the WF of SPASTS has its neg- 
ative distribution in phase space. This indicates that the 
WF of SPASTS always has the negative values at the 
phase space center a = (7 = 0), which is different from 
the case of single-photon-subtracted STS with a condi- 
tion Tic < sinh^ A [33| , but similar to the case of single- 
photon-added/subtracted squeezed vacuum [28l, [37j. 

Using Equations ([321)- ([Ml) we show the plots of WF 
in the phase space in Fig. 3 for the squeezing param- 
eter (A = 0.3) with different photon-added numbers m 
and average number Uc of the thermal state. One can 
see clearly that there is some negative region of the WF 
in the phase space which implies the nonclassicality of 
this state. In addition, the squeezing effect in one of the 
quadratures is clear in the plots (see Figure 3(a)), which 
is another evidence of the nonclassicality of this state. 
The WF has its minimum value for m = 1, 3 at the cen- 
ter of phase space (a — 0) (see Fig. 2(a) and (d)). The 
case is not true for m — 2 (see Fig. 2(c)). For m — 2, 
there are two negative regions of the WF, which differs 
from the case of single PASTS. In addition, the negative 
region of WF gradually decreases with the increasement 
of Uc, but not disappear for m = 1. 
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FIG. 3: (Color online) Wigner function distributions 
W {a, a*) of PASTS with A — 0.3 for different ric and m values 
(a) ric — 0.1, m = l;(b) ric — 0.5, m — 1; (c) Uc — 0.1, m — 2; 
(d) ric = 0.1, m = 3. 



V. DECOHERENCE OF PASTS IN THERMAL 
ENVIRONMENT 

In this section, we consider how this single-mode 
PASTS evolves at the presence of thermal environment. 
In thermal channel, the evolution of the density matrix 
for the m-PASV can be described by [11] 

— = K (Af + 1) (2apa^ — a) ap — pa) a) 
at 

+ kM (2a^ pa — aa^ p — paa^) , (38) 

where k represents the dissipative coefficient and Af de- 
notes the average thermal photon number of the envi- 
ronment. When Af = 0, Eq. (|38|) reduces to the master 
equation describing the photon-loss channel. The evolu- 
tion formula of WF of the PASV can be derived as follows 

M 

(39) 

where II^(Q;,a*,0) is the WF of the initial state, and 
T = 1 - e"^"*. Thus, in thermal channel, the WF at 
any time can be obtained by performing the integration 
when the initial WF is known. 

In a similar way to deriving Eq. p2p . substituting 
Eqs. ((32| - (p4| into Eg. ([39)) and using the generating func- 
tion of single- variable Hermite polynomials ([S|) , we finally 
obtain (see Appendix B) 
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Equation (PH]) is just the analytical expression of WF for 
PASTS in the thermal channel. It is obvious that the 
WF loses its Gaussian property due to the presence of 
single-variable Hermite polynomials. It is interesting to 
notice that Wq {r],r]* ,t) is actually the WF of squeezed 
thermal state in thermal channel corresponding to the 
case without photon addition (m = 0), Fq {r],r]*,t) = 1; 
while Fm {r],rj* ,t) is just the non-Gaussian contribution 
from photon-addition. The partial negativity of WF is 
fully determined by that of Fm {ij, ij* ,t). 

In particular, at the initial time {t — 0), noting 
{2J\f+l)rVG 2, gl/G 1, and Ai ^ 250, A2 ^ 



-(cosh^ A + nc)/{2nc + 1), 
,2 



as 



sinh2A/[4(2nc + l)j, x 

well as uj/(2iy/A^) 7 = [77*sinh2A - 2r7(cosh^ A + 
nc)]/{i[i2nc + 1) smh2X\^/^}, Eqs.dHJ and ^ just do 
reduce to Eas. ([551) and (IMl) . respectively, i.e., the 
WF of the PASTS. On the other hand, when nt ^■ 
00, noticing that T l,G 4/ (2nc + 1)^ , Ai ^■ 
2/(2A/'+l),w/(2iVA^) ^ 0,A2 ^ i (2nc + l)sinh2A, 



— >■ ric cosh 2A -f cosh A, as well as Hm (0) 
f^rn,2j, then Eq. pO)) becomes W {r],r]* ,00) 



and X 

l/[7r(27V-)- l)]exp[-2|77|V(2A/'+ 1)], a Gaussian distri- 
bution, which is independent of photon-addition number 
m and corresponds to the WF of thermal state with mean 
thermal photon number A/". This indicates that the sys- 
tem state reduces to a thermal state with mean photon 
number Af after an enough long time interaction with the 
environment. 

In addition, for the case of to = 1, corresponding to 
the case of SPASTS, Eq. (gO]) just becomes 



(45) 



It is obvious that when Fi{r],r]*,t) < 0, the WF of 
SPASTS in thermal channel has its negative distribution 
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FIG. 4: (Color online) Wigner function distributions 
W{a,a*) of PASTS with m = 1, = 0.3 for different 
Af, X and K,t values (a) Af = 0.2, A = 0.3, = 0.05;(b) 
A/" = 0.2, A = 0.3, Kt = 0.2; (c) A/" = 0.2, , A = 0.8, Kt = 0.05; 
(d) AA = 2, A = 0.3, = 0.05. 



in phase space. At the center of phase space rj = rj* — 0, 
the WF of SPASTS always has the negative values when 



X < 0, i.e 

250+536" 



(2-g3e'«*)/(2go + 53e- 



252) < (note 



- 2g2 > 0) leading to the following condition: 



1, 2A/'- 
Kt < ntc — - In ——r- 

2 2M - 



1 



(46) 



which is independent of the squeezing parameter A and 
the average photon number Uc of thermal state, there 
always exist negative region for WF in phase space and 
the WF of PASTS is always positive in the whole phase 
space when nt exceeds the threshold value ntc- Due to 
this and from Eq. (|46p , we can see how the thermal noise 
shortens the threshold value of the decay time. Compar- 
ing to the time threshold value of SPSSTS [s^ with the 
identical squeezed thermal state to that of SPASTS, 



t\,tn 



2 



2n^ 



1 



— sinh A 



27V + 1 cosh 2A + sinh^ A 



one can find a difference of e^''*'= — 6^"*'== 



2nc (ric + 1) 



{2N + 1) {uc cosh 2A + sinh^ A) 



(47) 



(48) 



which implies that the decoherence time of SPASTS is 
longer than that of SPSSTS. In this sense, the photon- 
addition Gaussian states present more robust contrast to 
decoherence than photon-subtraction ones. 

In Fig.3, the WFs of PASTS with m = 1 and = 0.3 
are depicted in phase space for several different Af, A 
and Kt values. It is easy to see that the negative region 
of WF gradually diminishes as the time nt increases (see 
Fig.3 (a) and (b)). In addition, the partial negativity of 
WF decreases gradually as Af (or A) increases for a given 
time (see Fig.3 (c) and (d)). The squeezing effect in one 
of the quadrature is shown in Fig. 4(c). For the case of 
large squeezing value A and small Uc and values, the 
single-PASTS becomes similar to a Schodinger cat state. 
The WF becomes Gaussian with the time evolution. 



VI. NON-GAUSSIANITY MEASURE FOR 
PASTS 

As well known, non-Gaussian operators (such as 
photon-adding/subtracting) can improve the nonclassi- 
cality and entanglement between Gaussian states [l^[l3| . 
One reason of such an enhancement is their amount of 
non-Gaussianity [40,|4l|. Recently, an experimentally ac- 
cessible criterion has been proposed to measure the de- 
gree based on the conditional entropy of the state with 
a Gaussian reference [4^]. Therefore, it is of interest to 
evaluate the degree of the resulting non-Gaussianity and 
assess this operation as a resource to obtain non-Gaussian 
states starting from Gaussian ones. Noting that the STS 
can be considered as a generalized Gaussian state, thus 
the fidelity between PASTS and STS may be seen as a 
non-Gaussianity measure. For this purpose, we define 
the fidelity by ^ 



T = tr {p,p) /tr [pI] 



(49) 



where ps and p are the STS (a generalized Gaussian 
state) and the PASTS, respectively. 

Noticing tr (p^) = l/(2nc + 1), and using the formula 
(CI), we finally obtain (see Appendix C) 



Ca,7. 



K2 
A 



1/2 



Prr 



nl {nc + 1)^ 
(2n, 



B/VA 
(50) 



sinh^ 2A 



1) 



where 

nc{nc + l) 

Ki = — — — cosh 2A, K2 

zric + 1 

(51) 

Eg. ([50)) is just the analytical expression for the fidelity 
between PASTS and STS. It is obvious that when m = 
(without photon- addition), = 1. Comparing to the 
fidelity J^, between PSSTS and STS (59) in Ref.lHi], one 
can clearly see that 



T 



12 P„, [Hl^ 



Cs 



P,n B 



(52) 



where Z — — (2nc -I- 1) sinh X, H ~ 7icCosh2A + 
sinh^ A. Eg. ([5^ implies that the ratio to fidelities is 
just that to the normalization factors. In particular, for 
TO = 1 (the case of SPASTS), Eq.§U\i reduces to 



T Uc cosh 2A + sinh A 
J^s rir cosh 2 A 



cosh^ A 



<1, 



(53) 



from which one can see that J^<J^s, i-e., the amount 
of non-Gaussianity for SPASTS is larger than that for 
SPSSTS. 

This point is made clear in Fig. 5, in which the fidelity 
J- between PASTS and STS as the function of squeez- 
ing parameter A for different photon-addition number 
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FIG. 5: (Color online) The fidelity T between PASTS 
(PSSTS) and STS as the function of squeezing parameter A 
for different photon-addition number m = 0, 1, 2, 3(nc = 0.2). 



m. As a comparision, the fidelity Ts between PSSTS 
and STS is also shown in Fig. 5, from which one can see 
that the fidelity decreases as the increment of photon- 
addition/subtraction number m, as expected. The fi- 
delity T increases monotonously with the augment of the 
squeezing parameter A. However, the case is not true for 
the fidelity Ts- For a given m value, the fidelity T is al- 
ways smaller than the fidelity Ts within the region shown 
in Fig. 5. In this sense, the amount of non-Gaussianity for 
PASTS is larger than that for PSSTS. 



VII. CONCLUSIONS 

In this paper, we investigate the nonclassical proper- 
ties and decoherence of single-mode PASTS when evolv- 
ing under a thermal environment. Based on the fact 
that squeezed number can be considered as an Hermite 
polynomial excitation squeezed vacuum, the normally or- 
dering form of PASTS is directly obtained, from which 
one can expediently calculate some quasi-distributions, 
such as Q-, P- and Wigner function; And the normal- 
ization factor of PASTS is analytically derived, which is 
just proved to be an m-order Legendre polynomial of the 
squeezing parameter r and average photon number ric 
of the thermal state, a remarkable result. Furthermore, 
for any photon-added number m-PASTS, the explicit ex- 
pression of WF is derived, which considered as a product 
of the WF of STS in thermal channel and a non-Gaussian 
distribution resulting from photon-addition. It is shown 
that the WF of SPASTS always has the negative values 
at the phase space center, which is different from the case 
of SPSSTS with a condition ric < sinh^ A. Then the ef- 
fects of decoherence to the nonclassicality of PASTS in 
the thermal channel is also demonstrated according to 
the compact expression for the WF. The threshold value 
of the decay time corresponding to the transition of the 
WF from partial negative to completely positive definite 
is obtained for SPASTS at the center of phase space. It 
is found that the WF has always negative value for all 
parameters r, ric if the decay time nt < ntc = ^ In 2'f!fXi ' 



a larger value than that of SPSSTS. 

A comparison between the nonclassicality and deco- 
herence of PASTS and PSSTS shows that the photon- 
addition non-Gaussian states present more robust con- 
trast to decoherence than photon-subtraction ones, which 
may be due to the amount of non-Gaussianity for 
SPASTS is larger that that for SPSSTS. On the other 
hand, in the limit of vanishing squeezing and Uc — 0, the 
PASTS reduces to a single-mode Fock state, remaining 
non-Gaussian, while the PSSTS becomes Gaussian, as it 
reduces to the single mode vacuum. Entanglement eval- 
uation investigation for photon-subtracted/added two- 
mode squeezed thermal state is a future problem. 
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Appendix A: Derivation of WF ([32]) for PASTS 

Substituting Eq. p^ into Eq. (PT|) and using the inte- 
gration formula (jl8p . we have 

p^(c,,a*) = Izir^e2|"l' / ^|/3|2-exp[-(l + i?) 

T1T2 J TT^ L 

c 



T1T2 



ds'^dt'" 



■ exp 



C 



+ {2a + s) P* - {2a* + t) P + - (/3*2 + p"^) 



= Wq {a,a*)Fm {a, a 
where we have set 
^1 



s=t=0 

(Al) 



Wo {a, a*) 



nTlT2 



■ exp 



A2 {a' +a*') -2A3\a\ 



(A2) 



Q2m 

F,n {a, a*) = C~J, — — exp 



+ {A2a* - A^a) t -I- {A2a - AiO*) s]^^^^^ , 

(A3) 



and 
Al 

A2 

A^ 



1 



(1 + BY - C2 (2n, + 1) 
2C sinh2A 



2 ' 



- 1 = 



(1 + BY - C2 2n, + 1 
2{B + l) ^ cosh2A 

{l + Bf -C^ 
2{B + 1) 

{1 + Bf-C^ 



2nr + 1 ' 



nc+ cosh A 
2nc + 1 



(A4) 
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Substituting Eq.(A3) into Eq.(A2) yields Eq.®, i.e., 
the WF of squeezed thermal state. 

Further expanding the exponential term st included in 
(A3) into sum series, and using the generating function 
of single- variable Hermite polynomials [13] , 



Qn 

Hn{x) = — exp (2a;t - 



(A5) 



which leads to 

Qn 



exp {At + Bt^ 



(^iVBy Hn A/{2iVB) 
-iVsY Hn \A/{-2iVB' 



(A6) 



thus we can see 



X exp 



2H\ ds"^dt"^' 
1=0 



^arn / j 



0.1 



s=t=0 
Q2m 



1=0 



2H\ a7'97*' 9s™at" 



X exp 

_ A. 
^ 2? 



s=t=0 



1=0 



where 7 — ^20!* — A^a, and 7 = 7/(1^/^2), i-e. 

a* sinh 2A — 2q; (cosh^ A + tic) 



iy/{2nc + l)sinh2A 
Then using the recurrence relation of _ff„(x). 



Eq.(A7) becomes 



2^n\ 



■Hn-l{x), 



(A7) 



(A8) 



(A9) 



2H\ 



Q2l _ Q21 



22m a™ 



Z=0 



(m!)" (-2A4M2)' 
;![(m-0!]^ 



(AlO) 



Substituting Eq.(A4) into Eq.(AlO) yields Eq.dM]). Thus 
we complete the derivation of WF Ea. (|32p by combing 
Eqs. (A2) and (AlO). 



Appendix B: Derivation of WF (gOD for PASTS 
in thermal channel 

Substituting Eqs. ([221)- O into Eq.dMl), we have 



M^(?7,77*,t) = 



TT (2nc + 1) 



92 



X exp 



4 ^ 



91ST 



a 



■ exp 



(250 + 53 



+ (g^ff + 92S - 2giT) a 

+ {93V + 92T - 2gis) a* + 52 (a^ + OL*'^)]s^r=o 

(Bl) 

where we have set 



cosh 2A 
50 = ^3 = —r, 91 



92 = M 



2nc + 1 ' 
sinh 2 A 
2nc + 1 ' 



A4 Tic+cosh A 
~2 ~ 2nc + 1 ' 

26-"* 



53 



(2AA+i)r' 



(B2) 



Further using the integration p^. Eq.(Bl) can be put 
into the form 

(77, 77* , t) = F,„ (7?, r;* , t) T^o (^7, ^* , t) , (B3) 

where (ry, ?7*,t) is defined in Eg. ipi . and 



F^(7y,r7*,i) = C-„i 



92 



■ exp [A2 (s^ + r^) 



(B4) 

here (A2,a;,x) are defined in Eq. (j44p . In a similar way 
to deriving Eq. (|32|) . we can further insert Eq. (B4) into 
Eq. (I42|). 

Appendix C: Derivation of fidelity ()50p between 
PASTS and STS 

The fidelity (tr {psp)) can be calculated as the overlap 
between the two WFs: 

{a,a*)Wp{a,a*) , (CI) 

where Wq {a, a*) is the WF of squeezed thermal state ps- 
Using Eq. ([5^ we may express Eq.(Cl) as 

tr(p,p) = 47r J F^{a,a*)W^ {a,a*)d^a. (C2) 

Then employing Eqs. ([311) and (A2),(A3) as weh as the 
integration formula (ITS)) , we can put Eq.(C2) into the 
following form: 



tr (psp) 



4C: 



{2nc + 1) ds^d 
exp 



■ exp 



91ST 



2nc + 1 9s™ 9r 



450 \af + 2g2 {a^ + a*^) 
{92T -2gis)a*]^^^ 
exp [KiST + Ko (s^ +t2)] 



+ (.92 s - ^gir) a + {g2T - 2gis) a*]^^^^^ 

c-1 



s=r=0 ■ 

(C3) 
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where Ki is defined in Ea. ([5T|) . and 

2nl + 2nc + l . , „, 
Ko = — ttt: 7^ sinh2A. 

4 (271c + 1) 

Similarly to deriving Eg. ([20|) . we have 

■ exp [Ko {k^ + t^) + Kikt] 



and K2 = Kl — AK"^ given in Eq. ([5T|) . which leads to 
Eq.pj). 



(C4) 



Qs^Qrr' 



k=t=0 



(C5) 
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